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We present a method to measure the effective field contribution to spin-transfer-induced interac- 
tions between the magnetic layers in a trilayer nanostructure, which enables spin-current effects to 
be distinguished from the usual charge-current-induced magnetic fields. This technique is demon- 
strated on submicron Co/Cu/Co nanopillars. The hysteresis loop of one of the magnetic layers in 
the trilayer is measured as a function of current while the direction of magnetization of the other 
layer is kept fixed, first in one direction and then in the opposite direction. These measurements 
show a current-dependent shift of the hysteresis loop which, based on the symmetry of the magnetic 
response, we associate with spin-transfer. The observed loop-shift with applied current at room 
temperature is reduced in measurements at 4.2 K. We interprete these results both in terms of a 
spin-current dependent effective activation barrier for magnetization reversal and a spin-current de- 
pendent effective magnetic field. From data at 4.2 K we estimate the magnitude of the spin-transfer 
induced effective field to be ~ 1.5 x 10~^ Oe cm'^/A, about a factor of 5 less than the spin-transfer 
torque. 



I. INTRODUCTION 

A spin-polarized current may interact with the mag- 
netic moment of a thin-film nanomagnet causing the re- 
versal of magnetization (switching) or other magnetic 
excitations such as spin waves. These spin-transfer ef- 
fects are of fundamental importance in understanding 
spin-transport at the nanoscale and are also of techno- 
logical relevance in the development of a new class of 
spin-electronic devices such as current-switched magnetic 
memories and current-tunable microwave sourcesi^i. 

The general form of the spin-transfer torque, intro- 
duced in 1996 by Slonczewski^ and Berger— , has found 
strong support in experiments^^SiLSiSiiSiiiiiS. Alternate 
models have been proposed, however, which consider an 
additional "effective field" interaction between the mag- 
netic moments in a multilayered structure. Heide et 
^^^13^14 proposed a model in which the effective field arises 
as a result of the longitudinal component of the spin ac- 
cumulation, which magnetically "couples" the layers to 
produce switching. An alternate diffusive model was pro- 
posed by Zhang et alM; in which the interaction between 
the local moments and the transverse component of the 
spin accumulation leads to the Slonczewski torque term, 
ajrh X (jh x nip), and an additional effective field term 
of the form bj m x nip, where rh and nip are unit vec- 
tors in the directions of the magnetization of the "free" 
and "fixed" magnetic layers, respectively. The coeffi- 
cients aj and bj both depend linearly on current and 
the ratio of their magnitudes depends on the thickness 
of the free layer and the decay length of the transverse 
component of the spin accumulation. In fact, any torque 
on the magnetization of the free layer can be decomposed 
into these two mutually perpendicular components in the 
plane orthogonal to rh. Brataas et alM- also considered 
a model with an effective field. In this model, the effec- 
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FIG. 1: Solid lines: Form of the zero-temperature stability 
diagram in the Slonczewski torque model for a single domain 
magnet with uniaxial anisotropy. Dashed lines: (a) Stability 
diagram with an additional effective field term of the form 
bjrhx nip, bj < 0. (b) Finite-temperature stability diagram 
in the Slonczewski torque model. In these schematic phase 
diagrams we have omitted the regions in which there are pre- 
cessional states for clarity—*—. 



tive field is related to the imaginary part of the mixing 
conductance G^^ while the spin torque is related to its 
real part. For metallic ferromagnets, ReG^-'^ is about ten 
times ImG'^^^, a similar conclusion to that reached by 
Slonczewski—. A related model by Waintal et ali^ and 
a model by Stiles and Zangwill et al?^ found no effective 
field type interaction. 

Within this framework, in the zero-temperature, mon- 
odomain approximation, the motion of the free layer 
magnetization under the infiuence of a spin-current is de- 
scribed by the Landau-Lifshitz-Gilbert (LLG) equation 
with the Slonczewski torque term and an effective field 
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term, 



dm 



Hcs - a 



dm 



(1) 



aj m X (to X mp) + bj m x mp. 



Hes is the effective field, = (iJa 



Hi,mx)x 



AirMm^z, where -ffappj the applied external field, and 
i?a, the in-plane uniaxial anisotropy field, are both in 
the X direction, z is the film normal and M is the free 
layer magnetization. 7 is the gyromagnetic ratio and a is 
the Gilbert damping parameter (a <C 1). The fixed layer 
magnetization top is also taken to be in the x direction. 
The coefficient in the Slonczewski term, aj, depends on 
the current density, the spin polarization P and the angle 
between the free and fixed magnetic layers rh ■ rhp. Fig- 
ure n^a) shows a schematic phase diagram obtained by 
numerically integrating Eq. ((SJ with the effective field 
term (dashed lines) and without it (solid lines). This 
effective field term leads to a linear current-dependent 
shift of the field axis, -ffapp ^^app — bj. The vertical 
boundaries which correspond to the external magnetic 
field overcoming the in-plane anisotropy field would 
be rotated clockwise or counterclockwise about J = 0, 
depending on the algebraic sign of bj. The spin torque 
term, however, does not effect this boundary at zero tem- 
perature. Below its threshold value, \aj\ < Gcq = a2-KM , 
the spin current modulates the damping^'^, but the mag- 
nitude of the damping does not change the switching field 
of the nanomagnet. Switching still occurs at the critical 
field when the metastable state becomes unstable. 

Finite temperatures may be modeled with the addi- 
tion of a Langevin random field Hl to the effective field 
H^s in Eq. This results in a finite probability 

for thermally activated switching with an activation bar- 
rier to magnetization reversal that depends linearly on 
the curren-fe2iiSS421. Figure ^fb) shows a schematic phase 
diagram where finite-temperature effects have been in- 
cluded. As can be seen from this figure and as we show 
below, at finite temperature the spin-torque term also 
leads to a shift of the midpoint of the hysteresis loop. 
However, if an effective field as that proposed by Zhang et 
al. were to be significant, its effects should be noticeable 
at low temperatures where the influence of thermal fluc- 
tuations becomes negligible. Experiments have ruled out 
the possibility of only a spin-transfer induced (STI) effec- 
tive field with no spin-transfer torque interaction, but not 
the possibility of both mechanisms concurrentlji2ii2L2&. 
Further, there have been no attempts to directly mea- 
sure the magnitude of the STI effective field interaction. 

It is important to note that charge-current induced 
(CCI) magnetic fields may also produce a shift of the 
hysteresis loop that is linear in the applied current. For 
example, if the contact resistance is comparable to that of 
the pillar, the current may flow asymmetrically produc- 
ing a net magnetic field bias on the nanomagnet. Thus, if 
we wish to measure a STI effective field, either one arising 
from a term of the form bj mxrhp or one associated with 
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FIG. 2: Differential resistance of a nanopillar spin valve device 
as a function of magnetic field measured at T = 295 K (a,b) 
and T = 4.2 K (c,d), for the fixed Co layer oriented in the 
+x (b,d) and —x (a,c) direction. The insets in (a,b) show the 
corresponding magnetic configurations of the layers, (c) Inset: 
Full-sweep differential resistance of the trilayer structure as a 
function of magnetic field measured at T = 4.2 K. Note that 
in this figure both layers eventually switch in the direction 
of the appfied magnetic field, while in (a,b,c,d) only the free 
layer switches its orientation. The vertical scale is the same 
as in (c,d). 



the Slonczewski torque term and finite temperature, we 
must identify the contribution from CCI magnetic fields. 

Here we describe a method to measure the effective 
field contribution of the spin-transfer mechanism that 
distinguishes the effect of a charge current from that of a 
spin current on the switching characteristics of the thin 
magnetic layer within a Co/Cu/Co trilayer nanopillar. 
Measurements at room temperature exhibit a linear shift 
of the center of the hysteresis loop of the thin magnetic 
layer with applied current, which is reduced in measure- 
ments at 4.2 K. We interprete these results in terms 
of both a Slonczewski torque term, including thermal 
fluctuations of the magnetization, and a spin-current- 
dependent effective magnetic field. Based on measure- 
ments at 4.2 K for which the influence of thermal fluctua- 
tions on the switching fields is negligible, we estimate the 
magnitude of a STI effective field of the form bjrhx rhp 
and compare this to the magnitude of the spin-torque 
term, aj m x (to x nip). 
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II. EXPERIMENTAL METHOD 

We study samples fabricated by thermal and elec- 
tron beam evaporation through a nano-stencil mask 
proces a^^'^" , with the stack sequence |3 nm Co|10 nm 
Cu|12 nm Co|300 nm Cu|10 nm Pt|. The thick "fixed" 
Co layer sets up a spin-polarization of the current, which 
then acts on the thin "free" Co layer causing the reversal 
of the orientation of its magnetization for large enough 
current densities (J w 10® A/cm^). Although the results 
we present below correspond to one sample of lateral size 
50 nm X 100 nm, similar results were obtained on sev- 
eral other samples with different lateral size dimensions. 
Transport measurements were conducted at 295 K and 
4.2 K in a four-point measurement geometry, where we 
measured the differential resistance dV/dl (and dc volt- 
age simultaneously) by means of a phase sensitive lock-in 
technique with a 200 fiA modulation current, at / = 750 
Hz, added to a dc bias current. The resistance per square 
of the top and bottom leads was found to be 0.65 fi, to 
be compared to the nanopillar resistance of 1.4 fi. The 
external magnetic field was swept at a rate of 12.8 Oe/s 
for measurements at T = 295 K and 20.3 Oe/s for mea- 
surements at T = 4.2 K. We define positive currents so 
that electrons flow from the fixed to the free Co layer. 

The orientation of the STI effective field is determined 
by the direction of the magnetization of the fixed Co 
layer. In particular, an effective field of the form bj m x 
rhp would tend to align the magnetization of the free 
layer m in the direction of that of the fixed layer rhp. 
Reversing the orientation of the fixed layer, i.e. mp 
^rhp, would simply change the sign of the effective field. 
On the other hand, CCI magnetic fields do not depend on 
whether the magnetic layers are aligned in the +x or —x 
direction. If Hg denotes the STI effective fields and He 
the CCI magnetic field, this symmetry may be written 
as 



Hs{mp) = BJ - 
Hc{mp) = CJ 



Hs{-rhp) = -BJ 
> Hs{-mp) = CJ, 



where Hg and He have both been expressed as linear 
functions of the current density. We thus see that in 
one magnetic configuration the CCI magnetic fields re- 
inforce the STI fields while in the other configuration 
they oppose the spin-transfer effects. Figure |21 shows 
these two magnetic configurations schematically together 
with the differential resistance versus magnetic field, with 
zero bias current, obtained for each configuration, at 295 
K (Figs. Ela) and (b)) and at 4.2 K (Figs. EJc) and 
(d)). Note the marked change in the width of the hys- 
teresis loop from 295 K to 4.2 K. The field drives the 
free layer moment hysteretically between a low-resistance 
state parallel (P) to the fixed Co layer and a higher- 
resistance antiparallel (AP) state. The magnetic cou- 
pling between the fixed and free Co layers (either ex- 
change or dipolar) causes the midpoint of the hysteresis 
loop to be shifted from H = 0. Carrying out magnetic 
field hysteresis loops for the free layer for both orienta- 
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FIG. 3: Phase diagrams of magnetic states for the free Co 
layer obtained from measurements of dV/dl as a function of 
H at fixed / (A, _ffp_Ap; V) Hap^p) and as a function of I 
at fixed H (O), measured at T = 295 K (a,b) and T = 4.2 K 
(c,d), for the fixed Co layer oriented in the +x (b,d) and —x 
(a,c) direction. Solid squares (■) show the observed effective 
(b,d) and HJ (a,c). 
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tions of the fixed layer, we are thus able to distinguish 
STI effective fields from CCI magnetic fields. The com- 
bined current-induced effective field on the nanomagnet 
can be decomposed into a charge-current and a spin- 
current component. By simply defining the observed to- 
tal effective field as 



-"P^AP 



(2) 



where the + or — corresponds to the effective field with 
the fixed layer in the +x or —x direction, it is easily 
seen from the preceeding discussion that the effective field 
contribution which only depends on the orientation of the 
fixed layer (STI) is given by 



Hg = 



H-, 



(3) 



and the contribution which is independent of the mag- 
netic configuration (CCI) by 



H-, 



(4) 



It is important to note that hysteresis loops must be per- 
formed on the free layer only. If we were to ramp the 
magnetic field so as to first switch the free layer and then 
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the fixed layer (Fig. EIc), inset), the effective field would 
be "reset" as we changed the orientation of the fixed layer 
and the free layer AP P field transition would not be 
measured. 



III. RESULTS 

In Fig. 13 we construct phase diagrams based on mea- 
surements of dV/dl as a function of H at fixed /, mea- 
sured at T = 295 K (Figs. Ela) and (b)) and at T = 4.2 
K (Figs. |2Ic) and (d)), for the fixed Co layer along the 
+x (Figs. Elb) and (d)) and -x (Figs. Eta) and (c)) di- 
rection. The figure shows at what values of H and / the 
different static states are stable or bistable. Note that, 
unlike previous measurementsSi, our phase diagrams do 
not include the switching of the fixed layer (See the inset 
in Fig. |2Ic)). Some of the switching boundaries are thus 
distinct from those obtained in field sweeps in which both 
fixed and free layers switch, in particular, the free layer 
AP — > P transition. Qualitatively, however, there are 
similarities between the two types of measurements, such 
as with respect to the effects of temperature. At T = 4.2 
K the field stability boundaries are almost independent 
of the applied current. The corresponding boundaries at 
T = 295 K, on the other hand, vary significantly with cur- 
rent. Also, while at T = 295 K the current dependence 
is clearly asymmetric, at T = 4.2 K the asymmetry is 
small. 

Figure 0] shows the different contributions to the to- 
tal effective field extracted from the data in Fig. 13 In 
Figs. Efb) and (c) we have used Eqs. l|2Jl and Q), re- 
spectively, after substracting the magnetic coupling field 
(determined from the shift from iJ = of the midpoint 
of the hysteresis loop for J = 0), to deduce the contri- 
butions from the fields that depend on magnetic config- 
uration (STI) and those that do not (CCI). Figure ^a) 
shows the current dependence of the width of the hys- 
teresis loop for T = 295 K and T = 4.2 K. 

From the data in Fig. 01b) it is clear that the effective 
field Hs is strongly dependent on temperature. While 
CCI magnetic fields are small both at room temperature 
and at 4.2 K (except for a clear effect for | J| > 0.3 x 10* 
A/cm^, which corresponds to the appearance of magnetic 
domains), the strong current dependence of the effective 
field at room temperature is reduced at 4.2 K. This data 
suggests that thermal fiuctuations play an important role 
in determining the switching boundaries. We thus con- 
sider the influence of finite temperature before estimat- 
ing the magnitude of the STI effective field interaction in 
these structures. 



IV. FINITE TEMPERATURE EFFECTS 

Finite temperature has been shown to have an im- 
portant effect on the spin-transfer induced switching 
thresholds2iiSSi^. In the low current regime or subcriti- 
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FIG. 4: Solid symbols denote measurements at T — 295 K 
and open symbols denote measurements at T = 4.2 K. (a) 
Hysteresis loop width as a function of current density for the 
fixed Co layer oriented in the +x (t>, ►) and —x (<1, m) direc- 
tion. The dashed line for T — 295 K was obtained by fitting 
the data near J = to Eq. ©, with A = [akBT/UoY^I^ as 
the fitting parameter. The dashed line for T = 4.2 K is a plot 
of Eq. Q with the value of Uo obtained from the T = 295 
K fit, /3 = 3/2 and a = 20. (b) STI effective field extracted 
from the data in Fig. lousing Eq. I^. The dashed lines are 
plots of Eq. O with the parameter Uo obtained from the fit 
to the T = 295 K hysteresis loop width data, j3 = 2i/2 and 
a = 20. The solid line is a linear fit to the STI effective field, 
near J = 0, at T = 4.2 K. (c) CCI magnetic field obtained 
from Fig. fusing Eq. @. 



cal region (J < Jco), where Jco is the zero-temperature 
threshold current density, finite temperature gives a non- 
vanishing probability for thermally activated switching. 
For any thermally activated process we may define an 
effective barrier U by fitting the mean switching rate to 
the form = exp{—U/kBT). The form of the ef- 
fective barrier has been shown to be approximated well 
by a power law, 



U, = Uo{lTKf{lTj), 



(5) 



where a — AP (upper sign) or P (lower sign) labels the ac- 
tivation barrier for transitions out of the AP or P states, 
Uq = (l/2)miJa, with m the magnetic moment of the 
nanomagnet, and ha = H^/ Hco = H„/ Hg. and j — J/ Jco 
are dimensionless variables corresponding to the exter- 
nal magnetic field and current, respectively, rescaled by 
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their zero-temperature critical values. The exponent 
/3 = 2 for the external magnetic field applied parallel 
to the easy axis and (3 = 3/2 for the external field ap- 
plied 10° away from the easy axis^. For a fixed value 
of j the switching field is determined by the attempt 
frequency and measurement time and thus depends on 
the applied magnetic field ramp-rate (r = 1 sec), i.e. 

w Tq-i exp{~U/kBT), which for t^^ = IQ-^ s-\ oc- 
curs when U/kBT ~ a k, 20. Using U — aksT for the 
barrier height we thus obtain the switching field 



h, = ±<l- 



akBT 



i/fr 



(6) 



Using Eq. Q we calculate the thermal-activation effec- 
tive field to be 



1 / aksT 

7 \W 



1//3 



3 



(l_j-2)l//3- 



(7) 



The width of the hysteresis loop, \hAp — hp\, as a function 
of current is given by 



^width 2 



1 - 



/ akBT\ 



1 



(8) 



neglecting higher order terms in j since we are interested 
in the regime where j <^ I. The finite-temperature field 
stability boundaries determined by Eq. © are shown 
schematically in Figure ^b). The trends are consistent 
with the phase diagrams shown in Figure |31 Figures 2Jb) 
and (c) show the results of fitting Eqs. {Tj) and (jS)), with 
J^o = 1.3 X 10® A/cm^ and Hco = = 320 Oe, to the 
bias field and the hysteresis loop width, respectively, in 
a way we detail next. In this treatment we are inter- 
ested in estimating the temperature contribution to the 
observed STI effective field. We thus use the hysteresis 
loop width data, which is unchanged by the presence of 
any effective field, to determine the energy barrier. Fit- 
ting Eq. |(SJ) to the hysteresis loop width as a function 
of current, at T = 295 K, in Fig. ^a), yields the value 
0.67 for the fitting parameter A = (ofesT/C/o)^/". If we 
take (3 = 3/2 and a = 20 this gives an energy barrier of 
0.93 eV. This is less than the value expected in a single 
domain picture of Uq ~ 2.2 eV, for M = 1440 emu/cm'^ 
and i?a = 320 Oe, but not unreasonable, as the mag- 
netization reversal likely occurs non-uniformly2^. Using 
these parameters, we have plotted Eq. ((SJ for the hys- 
teresis loop width at T = 4.2 K, in Fig. 2{a), and Eq. 
O for the effective field at T = 295 K and at T = 4.2 K, 
in Fig. Elb). The value we obtain for the energy barrier 
accounts for the scaling with temperature of the zero- 
current hysteresis loop width. Note that the measured 
width of the hysteresis loops decreases more rapidly than 
expected based on this thermal activation model, likely 
because the magnetization reversal processes are current 
dependent. For example, this current dependence could 
result from the circular field generated by the charge cur- 
rent influencing the element's domain structure. It is also 



clear from these curves that these fits underestimate the 
change in the effective field with current, both at T = 295 
K and at T = 4.2 K. At T = 4.2 K the effects of thermal 
ffuctuations are negligible near zero current, because the 
thermal energy is much less than the energy barrier to re- 
versal. This is evidence for an effective field of the form 
bj = m X rhp, in addition to the spin-transfer torque 
interaction. 



V. EFFECTIVE FIELD INTERACTION 

We can estimate the magnitude of a STI effective field 
of the form bj = rh x rhp based on the measurements 
at r = 4.2 K. If we account for the (relatively minor) 
contribution of thermal ffuctuations at T = 4.2 K to the 
current-depend shift in bias field, we obtain, near J = 
(fit in the range, -0.25 J^o < J < 0.25 J^o), -1.5 x 10^^ 
Oe cm^/A for the slope of bj. It is of interest to compare 
this to the magnitude of the spin-transfer torque inter- 
action. From Eq. Q, the spin-transfer-torque-driven 
switching boundary for the AP — > P transition is given 
by Oc = a(— i?app + H^. + 27rM), while that for the P 
AP transition is given by flc = — a(i?app + ^^a + 2ttM). 
If we take a = 0.01, we find 8.0 x lO"'^ Oe cm^/A, based 
on the measured switching current at zero applied field. 
To compare this quantitatively to the models^^'^'^ we 
compute the ratio of the magnitude of the spin-transfer 
torque aj to that of the effective field term bj. We find 
\aj/bj\ « 5.3. The values in the detailed diffusive trans- 
port model of Shpiro et alm^ are extremely sensitive to 
the Co layer thickness, the interface resistance and the 
decay length of the transverse component of the spin ac- 
cumulation, Aj. Within this model values of Xj ~ 2nm 
can produce a ratio close to that found in experiment. 
We note that the relative signs of these interactions are 
consistent with the models of Zhang et al. and Shpiro et 
al. 



VI. CONCLUSIONS 

In summary, we find evidence for a STI effective field 
interaction, with a magnitude that is a factor of 5 less 
than the spin-transfer torque for a 3 nm thick Co free 
layer. Measurements of the interaction constant bj as a 
function of free layer thickness will be essential to under- 
standing the origin of this interaction, such as whether 
it is associated with the longitudinal^"' or transversei^i^^ 
component of the spin accumulation. In the latter case, 
the interaction strength is expected to depend strongly 
on the free layer thickness, when this is comparable to 
the decay length of the transverse component of the spin 
accumulation, ^ 2 n"^-'^^ 



The longitudinal spin accu- 
mulation length is much longer (i.e. 60 nm in Co) and 
thus these effects should thus be easy to distinguish. Fu- 
ture studies will also examine a wider variety of device 
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structures and sample to sample fluctuations in these in- 
teraction constants. 

While a number experiments have been able to ex- 
clude the presence of a STI effective field without a spin- 
transfer torque iinteraction^-'^'^^-^* , these results show 
that both interactions are likely to be present. Note that 
for typical nanopillar samples and comparable interac- 
tion strengths, the current induced switching threshold is 
still mainly determined by the spin-transfer torque inter- 
action, while, as we have shown, at sufficiently low tem- 
perature and currents the magnetic switching threshold 
reflects the STI effective field interaction. The method 
we have presented enables a clear separation of these 



two mechanisms and can be applied quite generally to 
nanopillar devices with patterned fixed and free magnetic 
layers. 
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